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Abstract 

Using a non-perturbative method based on the selfconsistent Quasi-particle Random-Phase Ap- 
proximation (QRPA) we describe the properties of the pion in the hnear a-modeL It is found that the 
pion is massless in the chiral hmit, both at zero- and finite temperatm'e, in accordance with Goldstone's 
theorem. 



1 Introduction 



Chiral Lagrangians, as the effective low-energy realization of QCD, have become increasingly important 
in hadronic physics. In the sector of up and down quarks and for vanishing quark masses, QCD exhibits 
an exact global SU{2)l x SU{2)ji chiral symmetry. In the non-perturbative vacuum, this symmetry 
is spontaneously broken to its vectorial subgroup SU{2)y with the appearence of pions as Goldstone 
bosons. The interaction of the Goldstone bosons is greatly restricted by chiral symmetry involving the 
ratios (m,r/4'7r/^)^ and (ii^^/47r/7r)^ as a small expansion parameters, where m^r, E.,^ and f-j^ denote the 
pion mass, the pion energy and the weak pion decay constant respectively. A systematic expansion is 
provided by chiral perturbation theory For example, the vr — vr scattering amplitude is determined 
order by order in the number of derivatives. In this way the low-energy theorems are known to be 
maintained. For many reasons, however, it would be interesting to have a non-perturbative approach 
while still maintaining the low-energy theorems. One obvious reason is the requirement of unitarity of the 
S-matrix in the scattering problem. Another relates to the thermodynamics of effective chiral theories, 
especially in the study of chiral restauration. As the critical point is approached one cannot expect 
perturbation theory to provide a valid description. 

Needless to say that the question of preserving the symmetries non-perturbatively is a very delicate 
one . While in perturbative calculations the class of diagrams that ought to be considered in order to 
preserve the symmetries in the physical observables is well known, the situation is far less clear in the non- 
perturbative case. The aim of the present paper is to demonstrate that such a program is indeed possible. 
Our theoretical framework will be the linear cr-model which is especially suited for the techniques to be 
employed. These techniques have their origin in many-body physics and consist of a mean-field treatment 
via a Bogoliubov rotation supplemented by RPA fluctuations. It is well known that such an approach, 
while being non-perturbative, treats symmetries and spontaneous symmetry breaking correctly |^]. We 
shall demonstrate that, exactly as in the fermionic case, the RPA built on the selfconsistent mean field 
is able to restore the symmetry broken by the mean-field vacuum. 

The paper is organized as follows: First the formulation of the bosonic mean-field problem will be 
given in sect. 2. In the 'quasi-particle basis' thus obtained, the RPA excitation spectrum for the single- 
pion mode is constructed in sect. 3. It will be shown explicitly that this spectrum contains a zero mode 
in the chiral limit, to be identified with the 'Goldstone pion'. In sect. 4 the formalism will be extended 
to finite temperature as a first step towards a non-pertubative description of the chiral phase transition. 
Again, there is no mass generation in the chiral limit. Conclusions and an outlook are given in sect. 5. 
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2 The Bogoliubov Rotation 

The starting point is the Lagrangian density of the Unear cj-model 
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where Ao represents the bare coupUng constant, the mass parameter and tt and a denote the bare 
pion and sigma fields, respectively. Chiral symmetry is explicitly broken (in the PCAC sense) by the last 
term in the Lagrangian, ca. At tree level the pion and sigma masses are given by 



Mo + 3Ag(<7)2 



(2) 



The pion possesses manifestly the Goldstone boson character since its mass is trivially proportional to c. 
The perturbative one loop calculation preserves this result as is shown for instance in 
For further development it is now convenient to define the field operators in terms of creation and anni- 
hilation operators as 
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where the frequency ujq, common to both fields, is given by 



(4) 



In a first step a canonical transformation is performed for the pion as well as the sigma field. Thus 
we introduce a new set of creation and annihilation operators through the following Bogoliubov rotation 
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Uqdq Vqdq^ 



^qbq - Uqbq - Wq 



(5) 



with Uq, Vq, Xq and t/q being even functions of their argument, and Wq a c-number. The first equation is 
the usual bosonic Bogoliubov transformation applied to the pion field. The new vacuum with respect to 
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these Og operators is the well-known 'squeezed state'. In the second equation the transformation contains 
an additional 'shift' Wg to account for the macroscopic condensate {a). For later notation we will adopt 
the variable s to designate this condensate. To render the transformations canonical the Bogoliubov 
factors have to obey the following constraints 



ul-Vg = 1, 
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(6) 



In the 'quasi-boson' basis eq.(||) the fields and their conjugates read 
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and the quasi-boson vacuum |<1>) ( a|$) = f3\^) = ) is given by the following coherent state 



|<I>) = exp 



^zi{q)a+atg + Z2{q)b^btg + Wgbtg 



|0). 



(8) 



were |0) denotes the vacuum for the original basis (ag|0) = 6g|0) =0 ) and zi = ^, ^^2 = f • 

It is now straightforward to write the Hamiltonian of the linear sigma model in the quasi-particle basis. 

After normal-ordering one obtains 



H = noiv,y,s) + v[(3o+f3^] + [a+((?)aj(g)] + J^S.iq) [P+ {q)(3{q)] 



+ 



(9) 



where ":" in the interaction part of H denotes normal ordering (to avoid lengthy expressions the inter- 
action part is given in terms of field operators rather than in second-quantized form). 
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The pion and sigma fields are given by eq.(0), and the coefficients Ho, rj, <fo-,7r and Co-,7r read exphcitly 
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Here Jo and Jq are quadraticahy divergent integrals 
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arising from the tadpole loops in the selfenergies (see Fig. 1). 

As usual the amplitudes Ui , Vi and s are determined by minimizing the vacuum expectation value ^'^.^'^^ • 
This is in fact equivalent to demanding that the single-particle part of H be diagonal i.e. 07^,0- = 0, and 
that the term linear in the boson operators vanishes, i.e. rj = 0. Defined in this way the set of a and /3 
operators form the 'selfconsistent quasiparticle basis' (scqb). 

We now turn to the evaluation of the amplitudes u, v and x and y. First we note that the expressions 
for Ctt o- and £n a can be recast in the form 



c<s>iq) = {UqVq)e<s,{q) + 
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For notational purposes a generic field ^ has been introduced to designate either the pion or the sigma, 
and the corresponding Bogoliubov parameters (U,V) denote the pair (u,v) or (x,y). The following iden- 
tities are easily verified 



e^(g) = ujq + A^(g), 
Caiq) = LHq + Ao-(g), 
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With the above expressions and some trivial algebra one can extract selfconsistently the Bogoliubov 
factors from 



e$(g) - A$(g) 



4(g) -A|(g) 



(14) 



and the quasiparticle energies are given by 



(15) 



This result allows to reexpress the BCS gap equations for the auxiliary variables (U,V) in terms of more 
physical variables namely the quasi-pion and quasi-sigma masses as 



^^(0) 
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(16) 



In order to derive the BCS equations one should recall that we have made use of the two conditions 
Cn = Ca = arising from the minimization of TCoiv, y, s) with respect to v and y. The minimization with 
respect to s, on the other hand, yields an additional condition, namely r] = 0. This will fix the shift s via 

fil + Xy + 3Xl[Io + Jo] = ^. (17) 
The HFB results given above can be summarized diagrammatically as indicated in Fig. 1 



1 




Figure 1: Diagrammatic representation of the mean-field results for the pion and sigma. The dashed line denotes 
the selfconsistent quasi pion propagator, the solid line the quasi sigma, and the wavy line the two point Greens 
function of the bare field & of the Lagrangian density in eq. (|l|) . 

For a physical interpretation it is important to see now how the quasiparticle masses behave in the chiral 
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limit (c ^ 0). With the expressions given above the latter can be written as 

£l{^) = £ + 2A2[/o- Jo], 



s 

£l(Q) = -^+2Xls\ (18) 



which implies that the quasipion mass does not vanish due to the nonvanishing difference /q — Jo- This 
is in violation of Goldstone's theorem. To restore the symmetry one has to go further and this will be 
done in the next section. 

Before ending we wish to comment on the difference of the quasi-particle masses which can be written as 
and which arises from the finite value of the condensate, s, in the Goldstone phase. The explicit form of 



S7ro-(0) will be given later on. The expression eq.(19) is reminiscent of a Ward identity which links the 



three point function or tttto" vertex to the mass difference. An interesting feature of the identity above is 
that it only contains a weak divergence. By a simple redefinition of the coupling constant 

(Lo diverges only logarithmically) it can be rendered finite. It will turn out later that this redefinition 
of the bare coupling constant is also able to make the RPA solutions free of divergences. The full 
renormalization program will be discussed in a forthcoming paper. 

3 Single-Pion RPA 

Given the mean-field results presented in the last section the task is to obtain a Goldstone mode in the 
chiral limit. As is well known in many-body physics the restoration of a symmetry, which is broken at 
mean-field level, is provided by the "selfconsistent" RPA. To make this explicit for the case at hand we 
remind that Q'i\vac) represents the single-pion mode, where Q5 is the axial charge given by the volume in- 
tegral of the time component of the axial vector current. In the linear sigma model the current is given by 

Al^ = a{dV) - (a^a)7r^ (21) 
When expressed in the selfconsistent quasiparticle basis the axial charge then becomes 
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A remark is in order. We see that the operator Qf, when acting on the coherent state |$), as defined in 
eq. d^), can excite six different modes corresponding to a single-pion excitation and pairs of correlated 
pion and sigma excitations. When written in the original basis Q5 takes the same form as in eq. (|2 



except that the quasi-particle masses £-^,£(7 are replaced by the tree- level masses m.^,ma (eq. (|2D). It 
therefore leads to the same modes. The quasi-particle basis has the advantage, however, that in the chiral 
limit (c = 0) all modes survive, while in the original basis the single-pion excitation vanishes, since the 
tree-level pion mass goes to zero in that limit. We will see below that the quasiparticle representation is 
indeed needed. 

To proceed further, we consider the following RPA excitation operator 



(23) 



As usual the RPA ground-state correlations will be determined by the requirement that Qu\RPA) 
Applying the equation of motion method of Rowe [^, |^], one then has 



{RPA\ [5Qu , [H, Qt]] \RPA) = n^{RPA\ [6Q^ , Q+] \RPA) 



(24) 



which leads to the RPA equations for boons. In case of an exact symmetry, one particular solution has 
to be 'spurious', i.e. occurs at zero energy (Q^, = 0). To identify this solution one has to consider the 
operator which generates the symmetry. It has to be ensured, of course, that the latter possesses the 
excitations that are present in the general ansatz of the RPA operator Q'^ . Indeed one notices that the 
chiral symmetry operator Q^, when written in the original basis, has the same structure as the RPA 
operator. Two difficulties occur, however. The first has been eluded to and relates to the disappearance 
of the single-pion component from the symmetry operator when going to the chiral limit. The second 
difficulty is caused by the presence of the 'mixed' combinations ftgoj^"*" and b'^^atq in Qg. Such terms are 
undesirable since is no longer a solution of eq. (p^ . When written in the quasi-particle basis, the first 
problem is automatically cured, as mentioned above. The second, at first glance, seems to persist since 
the 'mixed' terms are still present. These terms give no contributions to the RPA equations , however, 
as long as the Hamiltonian is diagonal. By construction this is the case, of course. 

To make the spurious solution explicit, we consider the set of 4 coupled equations resulting from the 
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explicit form of the RPA excitation operator in eq. (23). Using Feshbach projection techniques it 
is advantageous to first solve the scattering problem for the pair of quasi-sigma and pion (lower part of 
Fig. 2) which is generated by the last two terms in . In the single-pion subspace one then has to solve 
a Dyson equation (upper part of Fig. 2) to finally obtain the physical pion mass. 



Figure 2: Upper part: The Dyson equation for the physical pion (thick dashed lines) for which the mass operator 
has been extracted from the scattering of the quasibosons in an RPA equation. Lower part: The scattering equation 
for a pair of quasi-sigma (thin full lines) and quasi- pion (thin dashed lines). 

This can in fact be done analytically and yields 

where 'Ena is the contribution of the (quasi) pion-sigma bubble to the pion selfenergy given by 

"""^ 7 (27r)3 2£^{q)£^iq) nl - {£^{q)+£^{q)y ^ > 



Using the identity 



we obtain after some algebra, 



o2 _ 2Ag [£^0) - £M] [S^.(O) - J^Unl)] c 

1 - 2XlE^Ani) + S- ^''^ 

In the chiral limit (c = 0) the zero-energy solution is now manifest (independent of any regularization 
scheme) . 

It is interesting to ask: What is the Goldstone boson dispersion relation? i.e. what is the behavior of the 
spurious mode under spatial translations. For this purpose we shall use as the generator of the spurious 
mode the time component of the axial vector current A^^{x) rather than the axial charge. The Fourier 
transform of the latter allows to pick up the spurious mode at any finite three momentum. To show that 
A^^(x) can generate a spurious mode we recall the PCAC relation 

a^^r(x)=C7r'^(x). (29) 
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Using Heisenberg's equation of motion, after Fourier transformation, this can be simply expressed as 



(30) 



In the chiral hmit the single-pion part of the RPA operator then generates a solution of finite three- 
momentum which has the following property 



{RPA\ 5Q, , H,Af{p) \RPA) + p{RPA\ 6Q, , A^^if) \RPA) = 0. 



(31) 



This clearly indicates that for pions at rest (p = 0) again a zero-frequency solution exist. To make the 
dispersion relation explicit we consider the following excitation operator 



q-p -q 



(32) 



which is the extension of eq. ( p^ ) to finite three-momenta. After some algebra the RPA frequencies can 
be expressed as 



2X^ [g^(0) - g^(0)] [S^.(O) - ^U^Ip)] 
1 - 2A2S^,(J^2^p) 
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with 



(2^)3 2£^{q)£^{p- q) - {£n{q) + £a{p - q)f 



(34) 



Using the fact that the vro" selfenergy is Lorentz invariant one easily verifies that 



^ilip) = nl{o) + f 

as it should be. Note that this result remains valid away from the chiral limit. 
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4 Finite Temperature HFB-RPA 

In a first step towards a nonperturbative description of the chiral phase transition we now extend the 
formalism to finite temperature by using well-known methods available in the literature|^, |9|, |l^. As 
we have demonstrated in the previous sections, the HFB-RPA has proven successful in preserving the 
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symmetry which is manifest through the presence of the spurious mode in the single pion RPA spectrum. 
The same is to be expected at finite temperature. Let us now first come to the mean field problem. We 
recall therefore that the thermodynamics of a gas of pions and sigmas at given temperature T is governed 
by the free energy Vt 

n = {H) - TS. (36) 

where (H) is the thermal expectation value of the Hamiltonian in eq. (^) and S denotes the entropy. In 
thermal equilibrium the distribution of maximum entropy is the one which minimizes fi. The entropy 
then reads 

S = kBY.[il + fu)Hl + /.) - fjnf,] (37) 
1/ 

where A:^ is the Boltzmann constant and fu is the usual bosonic distribution functions. The sum u 
includes the number of species as well as the three momentum q. 

In analogy to the zero-temperature case we can perform a temperature-dependent Bogoliubov rota- 
tion for both the vr and a fields. By making use of the Bloch-De Dominicis theorem |jll|] normal ordering 
on the rotated creation and annihilation operators can be carried out and one straightforwardly arrives 
at the HFB expression for the free energy 



HoivT, VT, St) = ^9(3^T,g + yT,g + 2) + 
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where 'Ho{vt , Ut , st) is just the expectation value of H on the grand canonical ensemble. Minimizing 
with respect to UT,q and XT,q and st while keeping the canonical normalization of the Bogoliubov factors 
as in eq. (^ then leads to the following identities 



r] 



SXqStIt + SXqStJt + AqSt + Mo-^t - c 
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where the definitions are as in the zero-temperature case and Ti.Q takes the same form as in eq. (|lO|). The 
loop integrals It and Jt are now given by 



It 



(27r)3 2uJq{uT,q - VT,qf 



Jt 
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1 + VM) 
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and the quasi-particle masses <Sj, <Sj and the condensate st take the form 



(40) 



^J(O)^ 
^J(O)^ 



5/r + Jt + 
It + Jt + 
fil + Ag4 + 3Ag (It + Jt) 



St 



(41) 



which is of identical form as the T = expressions. 

We now move on to the RPA problem at finite T. In the spirit of the zero-temperature RPA an oper- 
ator Q'^ is used which contains the same excitations as the symmetry generator. Therefore Q'^ is given by 

Q+=Xiag+ - + ^ [x2(g)/3+a"_+ - Y^\q)p.,a''^ + Xl{q)f3.,a-_+ - Y^{q)l3+af^ (42) 



which now contains additional terms of 'mixed' type. The analogous equations of motion have been 
worked out in refs. |l^] and read 



{[5Qu, [H,Qt]])=^A[^Qu,Qt]) 



(43) 



where the average is to be taken in the grand ensemble. The RPA operator (42) now generates a set of 
six equations for the various amplitudes which can be written as 



d^pMij{q,p)xj{p) = ^u^ij{q)xM) 



(44) 



where M.{q,p) is the 6x6 RPA matrix to be inverted, Vl^, and x{p) the eigenvalues and 6-column 
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eigenvectors respectively, and finally M{q) is the diagonal matrix norm 



ATuiq) = -N22{q) = l 

MM = -Mu{q) = I + Mq) + Mq) 

Mniq) = -Mmiq) = f'Al) - Mq) 



(45) 



The normalization condition for the eigenvectors is 



< 



+ iUq) - Uq)] (\ Xliq) P - I Y^{q) = 5,,, (46) 



The solution of the eigenvalue problem in eq.(|4J) gives 



(47) 



where 
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and the Bose occupation factors are given by 



fn,aip) 



expi 



£lAp). 



(49) 



One verifies that in the zero-temperature limit the previous HFB-RPA results are recovered. 
For the physical interpretation of these results we first address the question whether the FTHFB-RPA 
is able to describe the two realizations of the symmetry i.e. the Wigner phase and the Goldstone phase. 
In the Goldstone phase the vacuum condensate st is finite and the bare mass negative such that the 
third equation in (|4l| ) is satisfied in the chiral limit (c = 0). There exists therefore a Goldstone mode in 
the theory and the symmetry is no longer manifest in the particle spectrum which means that the masses 



12 



of the a and the vr are different. In the zero-temperature case we have demonstrated that HFB-RPA 
scheme is consistent with these requirements. 

There is, however, an alternative solution of the third equation in (|4l|). Suppose the condensate st 
vanishes at some temperature. In this case the three-particle coupling disappears from the interaction 
Hamiltonian, as can be seen explicitly from eq. (P). The single-particle state no longer couples to two- 
particle states and the only contribution of the single-particle masses is the one that comes from the 
four-point interactions in the mean-field calculation. This can also be checked explicitly from the RPA 
eigenvalues. It is now easy to see from eqs. (^Tj) that the masses become degenerate i.e. = S^^ and we 
are in the Wigner mode. The question is whether this is inconsistent with the fact that Q5, the generator 
of the symmetry, commutes with the Hamiltonian which leads to a spurious solution of the RPA in the 
Goldstone phase. First one should note that if the masses are equal then the thermal occupation factors 
for both the pion and the sigma are the same. This leads to A/55(g) = MeQ{q) = implying that the 
RPA spurious mode cannot be normalized. Secondly, from equation ( |47| ) we see that to have the zero 
frequency mode, one must fulfill the following condition 



In analogy to eq. (|27] ) one can prove the following identity 



which is only true, however, if all six terms in the excitation operator (eq. (|4^ ) are kept. Now the 
condition for a spurious mode solution can be simply recast as 

— = (52) 

ST 



which means that the ratio of the symmetry breaking term in the Lagrangian and the condensate must 
vanish to allow a zero-energy solution. This can only happen for finite st i.e. in the Goldstone phase. 
Once the Wigner phase is reached this condition can no longer be satisfied. This reiterates the fact that 
the spurious mode is a manifestation of a broken symmetry which disappears once the latter is restored. 
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5 Conclusions and outlook 



In summary we have presented a non-perturbative method based on the well-known selfconsistent QRPA 
formalism for studying the linear cj-model in the bosonic sector. Being 'symmetry conserving' the method 
yields a zero mode in the chiral limit. This is required by Goldstone's theorem for a spontaneously broken 
symmetry. While, at field level, the pion acquires a mass through the BCS mechanism irrespective of 
the explicit symmetry breaking term in the Lagrangian the inclusion of RPA correlations removes this 
artifact. We have also demonstrated that the extension of the QRPA to finite temperature is workable 
and reproduces the expected result that the zero mode persists at finite temperature. Applications to 
the 5C/(3)-case as well as the inclusion of fermions are straightforward and are being considered. This 
will hopefully provide new insight into the nature of the chiral phase transition. Since, in contrast to 
Nambu-Jona-Lasinio type models, the linear c-model is renormalizable a program of non-perturbative 
renormalisation should be persued in order to assess its impact on the physics. Despite the selfconsistency 
of the mean-field equations a solution of this problem does not seem out of reach. A challenging problem 
is the application of the formalism to the two-pion case. Here one attempts to build a scattering equation 
which is consistent with the low-energy theorems required by the symmetry. As is known from the 
analogous fermionic problem higher RPA schemes have to be employed. In particular the second RPA 
p^ , p^ ] is also 'symmerty conserving'. Its bosonic analog is easy to construct but some conceptual 
problems remain to be resolved. 
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